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In this work we discuss the place of Veneziano amplitudes (the precursor 
of string models) and their generalizations in the Regge theory of high energy 
physics scattering processes. We emphasize that mathematically such ampli- 
tudes and their extensions can be interpreted in terms of the Laplace (respec- 
tively, multiple Laplace) transform(s) of the generating function for the Ehrhart 
polynomial associated with some integral polytope V (specific for each scatter- 
ing process). Following works by Batyrev and Hibi to each such polytope V 
it is possible to associate another (mirror) polytope V. For this to happen, 
it is necessary to impose some conditions on V and, hence, on the generating 
function for V . Since each of these polytopes is in fact encodes some projective 
toric variety, this information is used for development of new symplectic and 
supersymmetric models reproducing the Veneziano and generalized Veneziano 
amplitudes. General ideas are illustrated on classical example of the pion-pion 
scattering for which the existing experimental data can be naturally explained 
with help of mirror symmetry arguments. 
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1 Introduction 



1.1 Brief history of the Veneziano amplitudes 

As is well known, the origins of modern string theory can be traced back to the 
4-particle scattering amplitude A(s,t,u) postulated by Veneziano in 1968 [1]. 
Up to a common constant factor, it is given by 

A(s, t, u) = V{s, t) + V(s, u) + V(t, u), (1) 

where 

l 

V(a,t) = J x-^-^l - x)- a ^dx = B(-a{s),-a{t)) (2) 
o 

is the Euler beta function and a(x) is the Regge trajectory usually written 
as a(x) = a(0) + a'x with a(0) and a' being the Regge slope and intercept, 
respectively. In the case of space-time metric with signature { — ,+,+,+} the 
Mandelstam variables s, t and u entering the Regge trajectory are defined by 
[2] 

s = -{pi+P2?] t = -(p 2 +p 3 ) 2 ; u=-(p 3 + Pl ) 2 . (3) 

The 4-momenta Pi are constrained by the energy-momentum conservation law 
leading to relation between the Mandelstam variables: 

4 

s + t + u = ^2m 2 . (4) 

i=l 

Veneziano [1] noticed 2 that to fit experimental data the Regge trajectories 
should obey the constraint 

a(s) +a(t) +a(u) = -1 (5) 

consistent with Eq.(4) in view of the definition of a(s). The Veneziano con- 
dition, Eq.(5), can be rewritten in a more general form. Indeed, let m,n,l 
be some integers such that a(s)m + a(t)n + a(u)l = 0. Then by adding this 
equation to Eq.(5) we obtain, a{s)m + a(t)h + a(u)l = —1, or more generally, 
a(s)rh + a(t)h + a(u)l + k ■ 1 = 0. Both equations have been studied extensively 

2 To get our Eq.(5) from Eq.(7) of Veneziano paper, it is sufficient to notice that his f — a(s) 
corresponds to ours -a(s). 
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in the book by Stanley [3] from the point of view of commutative algebra, poly- 
topes, toric varieties, invariants of finite groups, etc. Although this observation 
is entirely sufficient for restoration of the underlying physical model(s) repro- 
ducing these amplitudes, development of string-theoretic models reproducing 
such amplitudes proceeded historically quite differently In this work, we aban- 
don these more traditional approaches in favour of taking the full advantage of 
combinatorial ideas presented in Ref. [3] . This allows us to obtain models repro- 
ducing Veneziano amplitudes which are markedly different from those known 
in traditional string-theoretic literature. 

In 1967-a year before Veneziano's paper was published- the paper [4] by 
Chowla and Selberg appeared relating Euler's beta function to the periods of 
elliptic integrals. The result by Chowla and Selberg was generalized by Andre 
Weil whose two influential papers [5,6] brought into the picture the periods of 
Jacobians of the Abelian varieties, Hodge rings, etc. Being motivated by these 
papers, Benedict Gross wrote a paper [7] in which the beta function appears 
as period associated with the differential form " living" on the Jacobian of the 
Fermat curve. His results as well as those by Rohrlich (placed in the appendix 
to Gross paper) have been subsequently documented in the book by Lang [8]. 
Perhaps, because in the paper by Gross the multidimensional extension of beta 
function was considered only briefly, e.g. [7], p. 207, the computational details 
were not provided. These details can be found in our recently published papers, 
Rcfs. [9,10,11]. To obtain the multidimensional extension of beta function as pe- 
riod integral, following the logic of papers by Gross and Deligne [12], one needs 
to replace the Fermat curve by the Fermat hypersufrace, to embed it into the 
complex projective space, and to treat it as Kahler manifold. The differential 
forms living on such manifold are associated with the periods of Fermat hyper- 
surface. Physical considerations require this Kahler manifold to be of the Hodge 
type. In his lecture notes [12] Deligne noticed that the Hodge theory needs some 
essential changes (e.g. mixed Hodge structures, etc.) if the Hodge-Kahler mani- 
folds possess singularities. Such modifications may be needed upon development 
of our formalism. A monograph by Carlson et al, Ref. [13], contains an up to 
date exaustive information regarding such modifications, etc. Fortunately, to 
obtain the multiparticle Veneziano amplitudes these complications are not es- 
sential. In Ref. [10] we demonstrated that the period integrals living on Fermat 
hypersurfaces, when properly interpreted, provide the tachyon-free (Veneziano- 
likc) multiparticle amplitudes whose particle spectrum reproduces those known 
for both the open and closed bosonic strings. Naturally, the question arises: If 
this is so, then what kind of models are capable of reproducing such amplitudes 
? In this paper we would like to discuss some combinatorial properties of the 
Veneziano (and Vcncziano-like) amplitudes sufficient for reproducing at least 
two of such models: symplectic and supersymmetric. Mathematically, the re- 
sults presented below are in accord with those by Vergne [14] whose work does 
not contain practical applications. Before studying these models, we would like 
to make some comments about the place of Veneziano amplitudes and, hence, of 
whatever models associated with these amplitudes, within the Regge formalism 
developed for description of scattering processes in high energy physics. This is 
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accomplished in the next subsection. 



1.2 The Regge theory, theorem by Froissart, quantum 
gravity and the standard model 

As is well known, all information in particle physics is obtainable through proper 
interpretation of the scattering data.Thc optical theorem (see below) allows 
one to connect the imajinary part of the scattering amplitude with the total 
crossection a. By measuring this crossection experimentally one can obtain some 
information about the scattering amplitudes. Additional useful information 
can be obtained by collecting data for differential crossections, by using the 
dispersion relations, etc. [15]. There is an unproven common belief that in the 
limit of high energies all scattering processes are adequately described by the 
Regge theory [16, 17]. The Veneziano amplitude by design is Regge behaving 
[1]. To our knowledge, the proof that in the limit of high energies scattering 
amplitudes are Regge behaving had been obtained only for some special cases 
[16,17], including that of QCD [18]. Since, irrespective to their mathematical 
nature, all string theories are based on this ( generally unproven!) belief of the 
validity of the Regge theory, they can be as much trusted (even if totally correct 
mathematically !) as can be the Regge theory. 

In the Regge theory the experimental data are presented using the Chew- 
Frautchi (C-F) plot, Rcf.[16], pp. 144-145. On this plot one plots the Regge 
trajectories. Such trajectories relate particles with the same internal quantum 
numbers but with different spin (or angular momentum). From the standard 
string textbook, Ref.[2], it is known that for the open bosonic string the Regge 
trajectory is given by a(s) — a(0) + a's (in accord with Eq.(2) above). It is 
important though that a(0) = 1 and a' — 1/2 for the open string while a(0) = 2 
and a' = 1/4 for the closed string. In known string-theoretic formulations the 
numerical values of these parameters cannot be adjusted to fit the available 
experimental data since their values are deeply connected with the existing 
string-theoretic formalism [2] and, hence, are not readily adjustable. In the 
meantime, for high energies currently available it is known, e.g. read Rcf. [15], 
p. 41, that a(s) = 0.7 + 0.8s or a(s) = 0.44 + 0.92s for typical Regge 
tragectories. Claims made by some string theoreticians that the available range 
of high energies is not sufficient to test the predictions provided by the existing 
string theories cannot be justified because of the following. 

One of the major reasons for development of string theory, according to 
Rcf. [2] , lies in developing of consitent theory of quantum gravity. Indeed, in the 
case of closed bosonic string the massless (i.e. s = 0) spin two graviton occurs 
in the string spectrum only if a(Q) = 2. This fact alone fixes the value of the 
Regge intercept a(0) on the C-F plot to its value : a(0) = 2. As plausible as it 
is, such an identification creates some major problems. 

Indeed, in the case of 2 — > 2 scattering process the total cross section for the 
elastic scattering in s-channel (in view of the optical theorem, e.g. see Ref.[15], 
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p. 47) is given by 



cr(s) - s~ 1 ImA(s,t = 0), 



(6) 



where the scattering amplitude A(s, t) is either postulated (as in the case of 
Veneziano amplitude) or determined from some model (e.g. the standard string 
model [2], etc.). The above expression is valid rigorously at any energy. In the 
limit s — > oo the Regge theory provides the estimate for this exact result : 



where c is some constant. As it is with all processes described by the Regge the- 
ory [15-17], physically this result means the following : the analytical behaviour 
of the amplitude for elastic scattering in the s-channel is controlled (through 
the exponent in Eq.(7)) by the resonance in t— channel. In particular, if the 
resonance is caused by the graviton this leads the total crossection to behave as: 
<j{s) = cs. Unfortunately, the obtained result violates the theorem by Froissart. 
It can be stated as follows (e.g. see Ref.[16], p. 53) : 

Theorem 1.1. (Froissart) In the high energy limit : s — > oo the total 
crossection a(s) in s-channel is bounded by cr(s) s ^ 00 < const log(s/so) where 
so is some (prescribed) energy scale. 

Evidently, even if the current efforts (based on commonly accepted formal- 
ism) to construct mathematically meaningful string/brane theory eventually 
might succeed, such a theory will contradict the Froissart theorem for reasons 
just described. Hence, either this theorem is incorrect and should be recon- 
sidered or the underlying assumptions of string theory regarding gravitons are 
incorrect. 

Remark 1.2. The way out from this situation was recently developed in our 
recent work, Rcf. [18], where new equivalence principle for gravity is proposed 
based on known rigorous mathematical results. This new equivalence principle 
has major implications for the standard model of particle physics [19]. Since 
physical predictions based this model are in agreement with the Froissart the- 
orem already, the results of Ref.[18] effectively convert the existing standard 
model into a unified field theory accounting for all four types of known funda- 
mental interactions and being manifestly renormalizable and gauge- invariant. 

Incidentally, the intercept a(0) = 1 for the open string theory does have 
some physical significance. Indeed, in this case use of Eq.(7) produces a(s) = c' 
where d is yet another constant. Such high energy begaviour is typical for 
the pomeron-a, hypothetical particle like object predicted by Pomeranchuk- still 
undinscovered [15,17]. Additional ramifications of Pomeranchuk's work have 
lead to the prediction of the companion of the pomcron-the odderon [20] . 

In addition to the difficulty with the Froissart theorem, just described, the 
existing string-theoretic models suffer from several no less serious drawbacks. 



a(s) = cs' 



,a(t=0)-l 



(7) 
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For instance, the Regge theory in general and the Veneziano amplitude (a pre- 
cursor of the string model) in particular states that in addition to the lead- 
ing (parent) Regge trajectory there should be countable infinity of daughter 
trajectories-all lying below the parent trajectory. Nowhere in string-theoretic 
literature were we able to find a mention or an explanation of this fact. Ex- 
perimentally, however, typically for each parent trajectory there are only few 
daughter trajectories. In this work we shall provide a plausible theoretical ex- 
planation of this fact based on the mirror symmetry arguments. We would like 
to emphasize that since the models reproducing Veneziano amplitudes discussed 
below differ from those commonly discussed in string-theoretic literature, the 
numerical values for the slope a' and the intercept a(0) of the Regge trajectories 
can be readily adjusted to fit the experimental data. This is in accord with the 
original work by Veneziano [1] where no restrictions on the slope and intercept 
were imposed. 

1.3 Organization of the rest of this paper 

The rest of this work is organized as follows. Section 2 begins with some facts 
revealing the combinatorial nature of Veneziano amplitudes. This is achieved 
by connecting them with generating function for the Ehrhart polynomial whose 
properties are described in some detail in the same section. Such a polynomial 
counts the number of points inside the rational polytope (i.e. polytopc whose 
vertices are located at the nodes of the regular k— dimensional lattice) and at 
its boundaries (faces). In the present case the polytope is a regular simplex 
which is a deformation retract for the Fermat-type (hyper) surface living in 
the complex projective space [9,10]. Next, using general properties of generat- 
ing functions for the Ehrhart polynomials for the rational polytopes we discuss 
possible generalizations of the Veneziano amplitudes for polytopes other than a 
simplex. This allows us to use some results by Batyrev [22, 23] and Hibi [24] in 
order to introduce the mirror symmetry considerations enabling us to exclude 
the countable infinity of daugher trajectories on the C-F plot using mirror sym- 
metry arguments. General ideas are illustraded on the classical example of the 
pion-pion scattering [25] for which the existing experimental data can be natu- 
rally explained with help of mirror symmetry arguments. Next, in Section 3 we 
begin our reconstruction of the models reproducing Veneziano and the general- 
ized Veneziano amplitudes. It is facilitated by known connections between the 
polytopes and dynamical systems [14,26]. Development of these connections is 
proceeds through Sections 2-4 where we find the corresponding quantum me- 
chanical system whose ground state is degenerate with degeneracy factor being 
identified with the Ehrhart polynomial. The obtained final result is in accord 
with that earlier obtained by Vergne [14] whose work does not contain any phys- 
ical applications. In Section 5 the generating function for the Ehrhart polyno- 
mial is reinterpreted in terms of the Poincare' polynomial. Such a polynomial is 
used, for instance, in the theory of invariants of finite (pseudo)reflection groups 
[3,27]. Obtained indcntification reveals the topological and group-theoretic na- 
ture of the Veneziano amplitudes. To strengthen this point of view, we use some 
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results by Atiyah and Bott [28] inspired by earlier work by Witten [29] on su- 
persymmetric quantum mechanics. They allow us to think about the Veneziano 
amplitudes using the therminology of intersection theory [30] . This is consistent 
with earlier mentioned interpretation of these amplitudes in terms of periods of 
the Fermat (hyper)surface [9,10]. It also makes computation of these amplitudes 
analogous to those for the Witten -Kontsevich model [31, 32], whose refinements 
can be found in our earlier work, Ref.[33]. For the sake of space, in this work 
we do not develop these connections with the Witten- Kontsevich model any 
further. Interested reader may find them in Ref.[34]. Instead, we discuss the 
supersymmetric model associated with symplectic model described earlier and 
treat it with help of the Lcfshctz isomorphism theorem. This allows us to look 
at the problem of computation of the spectrum for such a model from the point 
of view of the theory of representations of the complex semisimplc Lie algebras. 
Using some results by Serre [35] and Ginzburg [36] we demonstrate that the 
ground state for such finite dimensional supersymmetric quantum mechanical 
model is degenerate with degeneracy factor coinciding with the Erhardt poly- 
nomial. This result is consistent with that obtained in Section 4 by different 
methods. 



2 The extended Veneziano amplitudes, the Ehrhart 
polynomial and mirror symmetry 

2.1 Combinatorics of the Veneziano amplitudes 

In view of Eq.(2), consider an identity taken from [37], 



= (1 + tz + (tz f + ...) •••(!+ tz n + (tz n f + ...) 



(1 - tz ) ■■■(!- tz k ) 



£( E *o*° ■■■*£*)*"■ (8) 

n— ko+... + k k — n 



When zq = ... = Zk = l,the inner sum in the last expression provides the total 
number of monomials of the type Zq" ■ ■ ■ with fco + ... + kk = n. The total 
number of such monomials is given by the binomial coefficient 3 

,. . (fc + n)! (n + l)(n + 2)---(n + fc) (fc + l)(fc + 2) • ■ ■ (fc + n) 
P(k,n) = = = : . 

(9) 

For this special case Eq.(8) is converted to a useful expansion, 

- oo 

P(k, t) = = ]T p(k, n)t n . (10) 

I 1 *) n=0 



3 The reason for displaying 3 different forms of the same combinatorial factor will be ex- 
plained shortly below. 
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In view of the integral representation of the beta function given by Eq.(2), we 
replace fc + 1 by a(s) + l in Eq.(lO) and use it in the beta function representation 
of the amplitude V(s,t). Straightforward calculation produces the following 
known in string theory result [2]: 

oo 

V(s,t) = -£p(a(s),n)- m —. (11) 

ra=0 ' 

The r.h.s. of Eq.(ll) is effectively the Laplace transform of the generating 
function, Eq.(lO). Such generating function can be interpreted as a partition 
function in the sence of statistical mechanics. 

The purpose of this work is to demonstrate that such an interpretation is 
not merely a conjecture and, in view of this, to find the statisical mechani- 
cal/quantum model whose partition function is given by Eq.(lO). 

Our arguments are not restricted to the 4-particle amplitude. Indeed, as we 
argued earlier [10,11], the multidimensional extension of Euler's beta function 
producing murtiparticle Veneziano amplitudes (upon symmetrization analogous 
to the 4-particle case) is given by the following integral attributed to Dirichlet 

V(x 1 ,...,x k )= J J n^- 1 nf-\..nf-\l-u 1 -...-u k ) x ^- 1 du 1 ...du k . 

ui>0,...,u k >0 
'Ui H \-u k <l 

(12) 

In this integral let t = u\ + ... + u k . This allows us to use already familiar 
expansion Eq.(10). In addition, the following identity 

t n = ( Ul + ... + u k ) n = V — X 1 (13) 

n=(n 1 ,...,n k ) 

with restriction n = n\ + ... + n k is of importance as well. This type of identity 
was used earlier in our work on Kontsevich-Witten model [33]. Moreover, from 
the same paper it follows that the above result can be presented as well in the 
alternative useful form: 



(«i + - + «k)" = f X Sx(ui, ...,« fc ), (14) 

Ahfc 

where the Schur polynomial S\ is defined by 

S A (Ul, .-.,«*)= Yl K ^n<---K k (15) 

n=(ni,...,n fc ) 

with coefficients K\ n known as Kostka numbers, / being the number of stan- 
dard Young tableaux of shape A and the notation Ahfc meaning that A is 
partition of k. Through such a connection with Schur polynomials one can de- 
velop connections with the Kadomtsev-Petviashvili (KP) hierarchy of nonlinear 
exactly integrable systems on one hand and with the theory of Schubert vari- 
eties on another. Although details can be found in our earlier publications 
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[33,11], in this work we shall discuss these issues a bit further in Section 5. Use 
of Eq.(13) in (12) produces, after performing the multiple Laplace transform, 
the following part of the multiparticle Veneziano amplitude 

A(l,...k)= . (a^+r)) 

(a(si) - m) • • • (a(sk) - n k ) 

Even though the residue F ni ...„ fe (a(sk+i)) contains all the combinatorial factors, 
the obtained result should still be symmetrized (in accord with the 4-particle 
case considered by Veneziano) in order to obtain the full murtiparticle Veneziano 
amplitude. Since in the above general multiparticle case the same expansion, 
Eq.(lO), was used, for the sake of space it is sufficient to focus on the 4-particle 
amplitude only. This task is reduced to further study of the expansion given by 
Eq.(lO). Such an expansion can be looked upon from several different angles. 
For instance, we have mentioned already that it can be interpreted as a partition 
function. In addition, it is the generating function for the Ehrhart polynomial. 
The combinatorial factor p(k, n) defined in Eq.(9) is the simplest example of the 
Ehrhart polynomial. Evidently, it can be written formally as 

p(k, n) = a n k n + a n ^ x k n - x + ■ ■ ■ + a . (17) 

2.2 Some facts about the Ehrhart polynomials 



A type of expansion given by Eq.(17) is typical for all Ehrhart-type polynomi- 
als. Indeed, let V be any convex rational polytope that is the polytope whose 
vertices are located at the nodes of some n— dimensional Z™ lattice. Then, the 
Ehrhart polynomial for the inflated polytope V (with coefficient of inflation 
k = 1,2,...) can be written as 

\kV nZ"|= $p(fc, n) = a n (V)k n + a^i^fc"" 1 + • • • + a (T) (18) 

with coefficients a , a„ being specific for a given type of polytope V. In the 
case of Veneziano amplitude the polynomial p(k, n) counts number of points in- 
side the n— dimensional inflated simplex (with inflation coefficient k = 1, 2, ...). 
Irrespective to the polytope type, it is known [38] that a = 1 and a n = VolV, 
where VolV is the Euclidean volume of the polytope. These facts can be eas- 
ily checked directly for p(k,n). To calculate the remaining coefficients of such 
polynomial explicitly for arbitrary convex rational polytope V is a difficult task 
in general. Such a task was accomplished only recently in [39] . The authors of 
[39] recognized that in order to obtain the remaining coefficients, it is useful to 
calculate the generating function for the Ehrhart polynomial. Long before the 
results of [39] were published, it was known [3,27], that the generating function 
for the Ehrhart polynomial of V can be written in the following universal form 

fe=o {i x) 



9 



The above result leading to possible generalizations/extensions of the Veneziano 
amplitudes does make physical sence as we shall demonstrate momentarily. Ad- 
ditional details are also presented in Section 5. 

The fact that the combinatorial factor p(k, n) in Eq.(9) can be formally writ- 
ten in several equivalent ways has physical significance. For instance, in particle 
physics literature, e.g. see [2], the third option is commonly used. Let us recall 
how this happens. One is looking for an expansion of the factor (1 — a;) - "^ -1 
under the integral of beta function, e.g. see Eq.(2). Looking at Eq.(19) one 
realizes that the Regge variable a(t) plays the role of dimensionality of Z- lat- 
tice. Hence, in view of Eq.(8), we have to identify it with n (or k, in case if 
Eq.(8) is used) in the second option provided by Eq.(9). This is not the way 
such an identification is done in physics literature where, in fact, the third op- 
tion in Eq.(9) is used with k = a(t) being effectively the inflation factor while 
n being effectively the dimensionality of the lattice 4 . A quick look at Eq.s(lO) 
and (19) shows that under such circumstances the generating function for the 
Ehrhart polynomial and that for the Veneziano amplitude are formally not the 
same. In the first case one is dealing with lattices of fixed dimensionality and 
is considering summation over various inflation factors at the same time. In 
the second case (used in physics literature [2]) one is dealing with the fixed 
inflation factor n = a(t) while summing over lattices of different dimensionali- 
ties. Nevertheless, such arguments are superficial in view of Eq.s(8) and (19) 
above. Using these equations it is clear that mathematically correct agreement 
between Eq.s(lO) and (19) can be reached only if one is using identification: 
<P(fc, n) — p(k, n), with the second option given by Eq.(7) selected. By doing so 
no changes in the pole locations for the Veneziano amplitude occur. Moreover, 
for a given pole the second and the third option in Eq.(9) produce exactly the 
same contributions into the residue thus making them physically indistinguish- 
able. The interpretation of the Veneziano amplitude as the Laplace transform 
of the Ehrhart polynomial generating function provides a very compelling rea- 
son for development of the alternative string-theoretic formalism. In addition, 
it allows us to think about possible generalizations of the Veneziano amplitude 
using generating functions for the Ehrhart polynomials for polytopes other than 
the n— dimensional inflated simplex used for the Veneziano amplitudes. As it 
is demonstrated by Stanley [3,27], Eq.(19) has a group invariant meaning as 
the Poincare' polynomial for the so called Stanley- Reisner polynomial ring. 5 . 
This fact alone makes generalization of the Veneziano amplitudes mathemat- 
ically plausible. From the same reference one can find connections of these 
results with toric varieties. In view of Ref.[14], this observation is sufficient 
for restoration of physical models reproducing the Veneziano and Veneziano-like 
generalized amplitudes. Thus, in the rest of this paper we shall discuss some 
approaches to the design of these models. 

Generalization of the Veneziano amplitudes is justified not only mathemat- 
ically. It is also needed physically as explained earlier in Subsection 1.2. The 

4 We have to warn our readers that nowhere in physics literature such combinatorial ter- 
minology is used to our knowledge. 

5 In Section 5 we provide some additional details on this topic. 
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information on Ehrhart polynomials just provided is sufficient for this purpose 
as we would like to explain now. 

2.3 The generalized Veneziano amplitudes and mirror sym- 
metry 

As we have explained already in Subsection 1.2., according to the Regge the- 
ory [16,17], for each parent trajectory there should be a countable infinity of 
daughter trajectories-all lying below the parent on the C-F plot. In his original 
paper [1], page 195, Veneziano took this fact into account and said explicitly 
that his amplitude is not uniquely defined. Following both the original work by 
Veneziano and Ref.[15], p. 100, we notice that beta function in Eq.(2) given by 
B{— a(s), —a(t)) (which is effectively the unsymmetrized Veneziano amplitude) 
can be replaced by B(m — a(s), n — a(t)) for any integers to, n > 0. To comply 
with the Regge theory one should use any linear combination of beta functions 
just described unless some additional assumptions are made. To our knowl- 
edge, the fact that the Veneziano amplitude is not uniquely defined regrettably 
is not mentioned in any of the existing modern string theory literature. Hence, 
if the alternative (to ours) formulations of string-theoretic models may finally 
produce some mathematically meanigful results, these formulations still will be 
confronted with explanation of the experimental fact that in nature only finite 
number of daughter trajectories is observed for each parent trajectory. If one 
accepts the viewpoint of this paper, such experimental fact can be explained 
quite naturally with help of mirror symmetry arguments. It should be noted, 
however, that our use of mirror symmetry differs drastically from that currently 
in use [40,41]. Nevertheless, the initial observations used in the present case 
do coincide with those used in more popular mirror symmetry treatments [41] 
since in our case they are also based on the work by Batyrev, Ref.[22]. In turn, 
Batyrev's results to some extent have been influenced by the result of Hibi [24] 
to be used in our work as well. 

Following these authors we would like to discuss properties of reflexive (polar 
(or dual)) polytopes. It is useful to notice at this point that the concept of the 
dual (polar) polytope was in use in solid state physics literature [42] for quite 
some time. Indeed both direct and reciprocal (dual) lattices are being used ruti- 
nely in calculations of physical properties of crystalline solids. The requirement 
that physical observables should remain the same irrespective to what lattice 
is used in calculations is completely natural. The same, evidently, should be 
true in the mirror symmetry calculations used in high energy physics. This is 
the physicall essence of mirror symmetry. In the paper by Greene and Plesser 
[43], p. 26, one finds the following statement : "Thus, we have demonstrated 
that two topologically distinct Calabi-Yau manifolds M and M' give raise to 
the same conformal field theory Furthermore, although our argument has been 
based only at one point in the respective moduli spaces Mm and Mm 1 of M 
amd M' (namely the point which has a minimal model interpretation and hence 
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respects the symmetries by which we have orbifolded) the results necessarily 
extends to all of Mm and Mm" ■ 

We would like to explan these statements now using more commonly known 
terminology. For this purpose we begin with the following 

Definition 2.1. A subset of R d is considered to be a polytope (or poly- 
hedron) V if there is a r x d matrix M (with r < d) and a vector b e R d 
such that P = {x£ R d | Mx < b}. Provided that the Euclidean c?-dimensional 

d 

scalar product is given by < x-y >= ^ x iV% > a rational ( respectively, integral) 

i=l 

polytope (or polyhedron) V is defined by the set 

P = {xeR d |<a,.x><ft ,i = l,...,r}, (20) 

where an g Q d (i-*) fc+1 and /3 i e Q for « = l,...,r (respectively ai e Z d and 
ft 6 Z for i = 1, ...,r\). 

Next, we need yet another definition 

Definition 2.2. For any convex polytope V the dual polytope V* is defined 
V* = {*E(R d ) | (a-x) < 1, ae V}. (21) 

Although in algebraic geometry of toric varieties the inequality (a • x) < 1 
is sometimes replaced by (a • x) > — 1 [38] we shall use the definition just stated 
to be in accord with Hibi [24]. According to this reference, if V is rational, then 
V* is also rational. However, V* is not necessarily integral even if V is integral. 
This result is of profound importance since the result, Eq.(19), is valid for the 
integral polytopes only. The question arises : under what conditions is the dual 
polytope V integral ? The answer is given by the following 

Theorem 2.3. (Hibi [24]) The dual polytope V* is integral if and only if 
T{T, x- 1 ) = {-l) d+1 xT(V, x) (22) 
where the generating function J-(V, x) is defined in Eq.(19). 

By combining Eq.s(10) and (19) we obtain the following result for the stan- 
dard Veneziano amplitude 

/ 1 \ d+1 

Wa0 =(_) . (23) 
Using this expression in Eq.(22) produces: 

T(V, x- 1 ) = = {-lf+^+^{V, x). (24) 

(1 - x) 
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This result idicatcs that scattering processes described by the standard Veneziano 
amplitudes do not involve any mirror symmetry since, as it is well known [22,23] 
in order for such a symmetry to take place the dual polytope V* must be in- 
tegral. In such a case both V and V* encode (define) the projective toric 
varieties Xp and Xpi which are mirrors of each other and are of Fano-type 
[22,23,45,46]. The question arises: can these amplitudes be modified with help 
of Eq.(19) so that the presence of mirror symmetry can be checked in nature? 
To answer this question, let us assume that, indeed, Eq.(19) can be used for 
such a modification. In this case we must require for the generating function 
T{V,x) in Eq.(19) to obey Eq.(22). Direct check of such an assumption leads 
to the desired result provided that h n -i = hi in Eq.(19). Fortunately, this is the 
case in view of the fact that these are the famous Dehn- Sommerville equations, 
Rcf.[38], p. 16. Hence, at this stage of our discussion, it looks like generalization 
of the Veneziano amplitudes which takes into account mirror symmetry is pos- 
sible from the mathematical standpoint. Unfortunately, in physics correctness 
of mathematical arguments is not sufficient for such generalization since exper- 
imental data may or may not support such rigorous mathematics. To check 
the correctness of our assupmtions (at least to a some extent) we would like 
to discuss now some known in literature results on pion-pion (ttit) scattering 
described, for example, in Refs. [25,46] from the point of view of results we just 
obtained. By doing so we shall provide the evidence that: a) mirror symmetry 
does exist in nature (wether or not its validity is nature's law or just a cu- 
riocity remains to be further checked by analysing the available experimental 
data) and, that b) use of mirror symmetry arguments permits us to eliminate 
the countable infinity of daughter trajectories allowed by the traditional Regge 
theory in favour of just several observed experimentally. 

Experimentally it is known that, below the threshold, that is below the 
collision energies producing more outgoing particles than incoming, the unsym- 
metrized amplitude A(s,t) for ttt: scattering can be written as 

A(s ' t] = - 92 T %-a(^-^)) )) = -9 2 (^(s)-a(t))B(l-a(s),l-a(t)). 

(25) 

This result should be understood as follows. Consider the "weighted" (still 
unsymmetrized) Veneziano amplitude of the type 

l 

A(s,t) = J dxx- a W- 1 (l-x)- a W- 1 g(x,8,t) (26) 
o 

where the weight function g(x,s,t) is given by 

g(x, s, t) = |.g 2 [(l - x)a(s) + xa(t)]. (27) 

Upon integration, one recovers Eq.(25). The same result can be achieved if , 
instead one uses the weight function of the type 

g(x,s,t) = g 2 xa(t). (28) 
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In early treatments of the dual resonance models (all developed around the 
Veneziano amplitude) [46] fitting to experimental data was achieved with some 
ad hoc prescriptions for the weight function g(x,s,t),e.g like those given by 
Eq.s (27) and (28). In the case of irir scattering such an ad hoc reasoning can 
be replaced by the requirements of mirror symmetry. Indeed, consider a special 
case of Eq.(19): n=2. For such a case we obtain, 

so that Eq.(22) holds indicating mirror symmetry. At this point, in view of 
Eq.s (26)-(29), one may notice that, actually, for this symmetry to take place in 
real world, one should replace the amplitude given by Eq.(25) by the following 
combination 

2 T(l-a( S ))T(l-a(t)) , 2 T(-a(s))T(-a(t)) 



A(s,t) = -g z \— w // 77^+3 



r(l - a{s) - a(t)) y T(-a{s) - a{t)) 
= -g 2 B(l-a(s),l-a(t))+g 2 B(-a(s),-a(t)). (30) 

Such a combination produces first two terms (with correct signs) of the infinite 
series as proposed by Mandclstam, Eq.(15) of Ref. [47]. The comparison with 
experiment displayed in Fig. 6. 2(a) of Ref. [46], p. 321, is quite satisfactory pro- 
ducing one parent and one daughter Regge trajectories. These are also displayed 
in Rcf.[15], p. 41, for the "rho family" of resonances. Thus, at least in the case 
of 7T7T scattering, one can claim that mirror symmetry consideration provides 
a plausible explanation of the observable data. One hopes, that the case just 
considered is typical so that mirror symmetry does play a role in Nature. 

The rest of this paper is devoted to the reconstruction of physical models 
reproducing Veneziano and extended Veneziano amplitudes based on mathe- 
matical results discussed in these two sections. Additional details can be found 
in Rcfs. [10, 11,34,48]. 



3 Motivating examples 



To facilitate our readers understanding, we would like to illustrate general prin- 
ciples using simple examples. We begin by considering a finite geometric pro- 
gression of the type 



T(c,m) = ^ expjd} = exp{— cm} ^ expje/} + exp{cm} exp{d} 



1=0 



1= — oo 



cxpl-cm} ^^ + expW ^^ ^ 
"exp{c(2m + 1)}-1" 



exp{— cm} 



exp{c} — 1 



(31) 
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The reason for displaying the intermediate steps will be explained shortly below. 
First, however, we would like to consider the limit : c — ► + of !F{c,m). It is 
given by ^"(0, m) = 2m + 1. The number 2m + 1 equals to the number of integer 
points in the segment [— m, m] including boundary points. It is convenient to 
rewrite the above result in terms of x = exp{c} so that we shall write formally 
T[x, m) instead of T[c, m) from now on. Using such notation, consider a related 
function 

fix, m) = (-l)^(-, -to). (32) 

x 

This type of relation (the Ehrhart-Macdonald reciprocity law) is characteris- 
tic for the Ehrhart polynomial for rational polytopes discussed earlier. In the 
present case we obtain, 

-(-2m+l) _ 1 

f(x,m) = (-l)- ; x m . (33) 

In the limit x — > 1 + + we obtain : f{l, m) = 2m — 1. The number 2m — 1 is 
equal to the number of integer points strictly inside the segment [— to, to]. Both 
J-{0,m) and J"(l,m) provide the simplest possible examples of the Ehrhart 
polynomials if we identify to with the inflation factor k. 

These, seemingly trivial, results can be broadly generalized. First, we replace 
£ by x —x\ ■ ■ ■ Xd, next we replace the summation sign in the left hand side of 
Eq.(31) by the multiple summation, etc. Thus obtained function J-(x, to) in the 
limit Xi — ► l+0 + , i = 1—d, produces the anticipated result : J-(l,m) = (2TO+l) d 
. It describes the number of points inside and at the faces of a d~ dimensional 
cube in the Euclidean space R d . Accordingly, for the number of points strictly 
inside the cube we obtain : f(l, to) = (2m — l) d . 

Let VertV denote the vertex set of the rational polytope. In the case consid- 
ered thus far it is the d— dimensional cube. Let {u\, u v d } denote the orthogonal 
basis (not necessarily of unit length) made of the highest weight vectors of the 
Weyl-Coxeter reflection group Bd appropriate for cubic symmetry [11]. These 
vectors are oriented along the positive semi axes with respect to the center of 
symmetry of the cube. When parallel translated to the edges ending at particu- 
lar hypercube vertex v, they can point either in or out of this vertex. Then, the 
c?-dimensional version of Eq.(31) can be rewritten in notations just introduced 
as follows 

cxp{< c • x >} = cxp{< c • v >} 

xG7>nZ d vEVertV 

(34) 

The correctness of this equation can be readily checked by considering special 
cases of a segment, square, cube, etc. The result, Eq.(34), obtained for the 
polytope of cubic symmetry can be extended to the arbitrary convex centrally 
symmetric polytope. Details can be found in Rcf. [49]. Moreover, the require- 
ment of central symmetry can be relaxed to the requirement of the convexity 
of V only. In such general form the relation given by Eq.(34) was obtained by 
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d 

n(l-exp{- Ci <}) 

_i=l 



Brion [50] . It is of central importance for the purposes of this work: the limiting 
procedure c — > + produces the number of points inside (and at the boundaries) 
of the polyhedron V in the l.h.s. of Eq.(34) and, if the polyhedron is rational and 
inflated, this procedure produces the Ehrhart polynomial. Actual computations 
are done with help of the r.h.s. of Eq.(34) as will be demonstrated below. 



4 The Duistermaat-Heckman formula and the 
Khovanskii-Pukhlikov correspondence 

Since the description of the Duistermaat-Heckman (D-H) formula can be found 
in many places, we would like to be brief in discussing it now in connection 
with earlier obtained results. Let M = M 2n be a compact symplectic manifold 
equipped with the moment map $ : M — > R and the (Liouville) volume form 
dV = (^r) 0™. According to the Darboux theorem, locally O = Ym=i dqi 
Adpi . We expect that such a manifold has isolated fixed points p belonging 
to the fixed point set V associated with the isotropy subgroup of the group G 
acting on M . Then, in its most general form, the D-H formula can be written 
as [14,26,51] 

dVe =Y,w^r> (35) 

where a\. pi a n , p are the weights of the linearized action of G on T p M. Using 
Morse theory, Atiyah [52] and others, e.g. see Ref.[54] for additional references, 
have demonstrated that it is sufficient to keep terms up to quadratic in the 
expansion of $ around given p. In such a case the moment map can be associated 
with the Hamiltonian for the finite set of harmonic oscillators. In the properly 
chosen system of units the coefficients a\ iP , a n>p are just "masses" rrii of 
the individual oscillators. Unlike truly physical masses, some of m-s can be 
negative. 

Based on the information just provided, we would like to be more specific 
now. To this purpose, following Vergne [53] and Brion [50], we would like to 
consider the D-H integral of the form 



I(k ; 2/1,2/2)= / dxidx 2 exp{-(yiXi + 2/2^2)}, 

JkA 



(36) 



where kA is the standard dilated simplex with dilation coefficient fc 6 . Calcula- 
tion of this integral can be done exactly with the result: 

I e~ kyi e~ km 

I(k ; y u 2/2) = + — r + — r (37) 

2/12/2 2/1(2/1 - 2/2) 2/2(2/2 - 2/1) 



6 Our choice of the simplex as the domain of integration is caused by our earlier made 
observation [10] that the deformation retract of the Format (hyper)surface (on which the 
Veneziano amplitude lives ) is just the standard simplex. Since such Fermat surface is a 
complex Kahlcr-Hodge type manifold and since all Kahler manifolds are symplectic [26,54], 
our choice makes sense. 
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consistent with Eq.(35). In the limit: 1/1,2/2 — * some calculation produces the 
anticipated result : VolkA = fc 2 /2! for the Euclidean volume of the dilated 
simplex. Next, to make a connection with the previous section, in particular, 
with Eq.(34), consider the following sum 

s(k ; 2/1,3/2) = E cx p{-(yi l i + y^)} 

{h,l 2 )ekA 

11 1 e-^ 1 1 e~ ky2 

~ 1 - e-f 1 1 - e-^ + l-e»l- e v^-yi \- e wi\ - e V2-yi 

(38) 

related to the D-H integral, Eq.s(36,37). Its calculation will be explained mo- 
mentarily. In spite of the connection with the D-H integral, the limiting proce- 
dure: 2/1,2/2 — ► in the last case is much harder to perform. It is facilitated by 
use of the following expansion 

1 1 1 s „ , o „ , 

7--7 = -+2 + 12- + ^- < 39 > 

Rather lcnghty calculations produce the anticipated result : S(k ; 0,0) = fc 2 /2! + 
3fc/2+l = |fcA n Z 2 \ = *P(fc, 2) for the Ehrhart polynomial. Since generalization 
of the obtained results to simplicies of higher dimensions is straightforward, the 
relevance of these results to the Veneziano amplitude should be evident. To 
make it more explicit we have to make several steps still. First, we would like 
to explain how the result, Eq.(38), was obtained. By doing so we shall gain 
some additional physical information. Second, we would like to explain in some 
detail the connection between the integral, Eq.(37), and the sum, Eq.(38). Such 
a connection is made with help of the Khovanskii-Pukhlikov correspondence. 

We begin with calculations of the sum, Eq.(38). To do this we need a 
definition of the convex rational polyhedral cone a. It is given by 

a = Z> ai H h Z> Q a d , (40) 

where the set ai,...,a<i forms a basis (not nesessarily orthogonal) of the d- 
dimcnsional vector space V, while Z> are non negative integers. It is known 
that all combinatorial information about the polytope V is encoded in the 
complete fan made of cones whose apexes all having the same origin in common. 
Details can be found in literature [26,30]. At the same time, the vertices of V 
are also the apexes of the respective cones. Following Brion[50], this fact allows 
us to write the l.h.s. of Eq.(34) as 

f(v,x)= E xm = E x<T ( 41 ) 

m£VnZ d aeVertV 

so that for the dilated polytope the above statement reads as follows [50,55]: 

n 

f(kV,x)= E x m = Ex kv *E x<Tl - ( 42 ) 

m£kVnZ d i=l c» 
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In the last formula the summation is taking place over all vertices whose location 
is given by the vectors from the set {vi,...,v„}. This means that in actual 
calculations one can first calculate the contributions coming from the cones 
<7j of the undilated (original) polytope V and only then one can use the last 
equation in order to get the result for the dilated polytope. 

Let us apply these general results to our specific problem of computation of 
S(k ; 2/1,2/2) in Eq.(38). We have our simplex with vertices in x-y plane given 
by the vector set { vi=(0,0), v 2 =(l,0), V3=(0,l)}, where we have written the 
x coordinate first. In this case we have 3 cones: u\ = Z 2 v 2 + Z 3 v 3 , er 2 = 
v 2 + Ji(-v 2 ) + k(v 3 -v 2 );<7 3 =v 3 + Z 3 (v 2 -v 3 ) + Zi(-v 3 );{Zi, l 2 , h }e Z> 
. In writing these expressions for the cones we have taken into account that, 
according to Brion, when making calculations the apex of each cone should be 
chosen as the origin of the coordinate system. Calculation of contributions to 
the generating function coming from o\ is the most straightforward. Indeed, 
in this case we have x = x\Xi — e~ vi e~ y2 . Now, the symbol x CT in Eq.(41) 
should be understood as follows. Since Oi , i = 1 — 3, is actually a vector, it has 
components, like those for vi, etc. We shall write therefore x CT = x^ 1 ^ • • • x a d ^ 
where a(i) is the i-th component of such a vector. Under these conditions 
calculation of the contributions from the first cone with the apex located at 
(0,0) is completely straightforward and is given by 

Y 4 2 4 3 = — - — — - — ■ (43) 

It is reduced to the computation of the infinite geometric progression. But physi- 
cally, the above result can be looked upon as a product of two partition functions 
for two harmonic oscillators whose ground state energy was discarded. By do- 
ing the rest of calculations in the way just described we reobtain S(k ; j/i,2/ 2 ) 
from Eq.(38) as required. This time, however, we know that the obtained result 
is associated with the assembly of harmonic oscillators of frequencies ±2/1 and 
±2/2 and ±(j/i — 2/2) whose ground state energy is properly adjusted. The "fre- 
quencies" (or masses) of these oscillators are coming from the Morse-theoretic 
considerations for the moment maps associated with the critical points of sym- 
plcctic manifolds as explained in the paper by Atiyah [52]. These masses enter 
into the "classical" D-H formula, Eq.s(36),(37). It is just a classical partition 
function for a system of such described harmonic oscillators living in the phase 
space containing critical points. The D-H classical partition function, Eq.(37), 
has its quantum analog, Eq.(38), just described. The ground state for such 
a quantum system is degenerate with the degeneracy being described by the 
Ehrhart polynomial ^3(fc, 2). Such a conclusion is in formal accord with results 
of Vergne [14]. 

Since (by definition) the coefficient of dilation k—1,2,... , there is no dynami- 
cal system (and its quanum analog) for k=0. But this condition is the condition 
for existence of the tachyon pole in the Vcncziano amplitude, Eq.(2). Hence, in 
view of the results just described this pole should be considered as unphysical 
and discarded. Such arguments are independent of the analysis made in Ref.[10] 
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where the unphysical tachyons are removed with help of the properly adjusted 
phase factors. Clearly, such factors can be reinstated in the present case as well 
since their existence is caused by the requirements of the torus action invariance 
of the Veneziano-like amplitudes as explained in [10,26]. Hence, their presence 
is consistent with results just presented. 

Now we are ready to discuss the Khovanskii-Pukhlikov correspondence. It 
can be understood based on the following generic example taken from Ref. [56] . 
We would like to compare the integral 

* f 

/e tz e sz v— ^ e tz e sz 

dxe zx = with the sum S(z) = > e kz = — h- . 
z z v ' ^ 1 - e- z 1 - e z 
s k=s 

To do so, following Refs[56-58] we introduce the Todd operator (transform) via 

Td{z) = r=W (44) 

Then, it can be demonstrated that 

f) f) ft+h 2 * 

rd tew^y._ hl e "^i^=°=g e **- (45) 

This result can be now broadly generalized. Following Khovanskii and Pukhlikov 
[57], we notice that 



Td (^) ex P [YlP iZ ij = Td(p!,...,p n )exp \^>2PiZiJ ■ (46) 



By applying this transform to 



1 

i{x 1 ,...,x k ;^ 1 ,.. .,4) = - — — exp(Va;i^) (47) 



we obtain, 



5(2;!,...,^;^,...,^) = — exp(^x^ 4 ). (48) 

IKl-exp(-O) 1=1 

i=l 

This result should be compared now with the individual terms on the r.h.s. of 
Eq.(34) on one hand and with the individual terms on the r.h.s of Eq.(35) on 
another. Evidently, with help of the Todd transform the exact "classical" results 
for the D-H integral are transformed into the " quantum" results of the Brion's 
identity, Eq.(34), which is actually equivalent to the Weyl character formula 
[48]. 

We would like to illustrate these general observations by comparing the D- 
H result, Eq.(37), with the Weyl character formula result, Eq.(38). To this 
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purpose we need to use already known data for the cones <7j , i = 1 — 3, and the 
convention for the symbol x CT . In particular, for the first cone we have already 
: x* 71 = x l ^x l 2 = [exp(7iyi)] • [exp(? 2 2/2)] 7 ■ Now we assemble the contribution 
from the first vertex using Eq.(37). We obtain, [exp(7iyi)] • [exp(/ 2 j/2)] /yiV2- 
Using the Todd transform we obtain, 

^(^)Td(^)^[exp to )].[exp( iaItt )] | w = 1 _[_ yi 1 ■ (49) 

Analogously, for the second cone we obtain: xf = e~ kyie~ l\y\e~ h(yi — 2/2) so 
that use of the Todd transform produces: 

Td{^-)Td{-^-)—^ e-ky^-hv^-Hv.-v,) \ h=h=a . ' 



y dh' y dh yi(yi-V2) m-n-u x _ eVl x _ eVl - V2 > 

(50) 

etc. 

The obtained results can now be broadly generalized. To this purpose we 
can formally rewrite the partition function, Eq.(24), in the following symbolic 
form 

I(k,f) = Jdxcxp(-f-x) (51) 

feA 

valid for any finite dimension d. Since we have performed all calculations ex- 
plicitly for two dimensional case, for the sake of space, we only provide the idea 
behind such type of calculation 8 . In particular, using Eq.(37) we can rewrite 
this integral formally as follows 



/ 



rf xcxp(-f.x)^ CXP( 7 f ' X(P)) - (^) 

i 

Applying the Todd operator (transform) to both sides of this formal expres- 
sion and taking into account Eq.s(49),(50) (providing assurance that such an 
operation indeed is legitimate and makes sense) we obtain, 



n(l-exp{-^(f)i 



,. d 

J /^n i _ ex p ( _ X4) -P(-f-)= E exp{<f.v>} 

fcA i=l veVertV 

= ]T cxp{<f-x>}, (53) 

where the last line was written in view of Eq.(34). From here, in the limit 
: f = we reobtain p(k,n) defined in Eq.(10). Thus, using classical partition 

7 To obtain correct results we needed to change signs in front of £1 and I2 . The same should 
be done for other cones as well. 

8 Mathcmatically inclined reader is encoraged to read paper by Brion and Vergne, Ref. [59], 
where all missing details are scrupulously presented. 
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function, Eq.(51), (discussed in the form of Exercises 2.27 and 2.28 in the book, 
Rcf.[58], by Guillemin) and applying to it the Todd transform we recover the 
quantum mechanical partition function whose ground state provides us with the 
combinatorial factor p(k, n). 

5 From analysis to synthesis 



5.1 The Poincare' polynomial 

The results discussed earlier are obtained for some fixed dilation factor k. In 
view of Eq.(8), they can be rewritten in the form valid for any dilation factor 
k. To this purpose it is convenient to rewrite Eq.(8) in the following equivalent 
form: 

det(l-Mt) = {l-t Z0 ) l ..(l-t Zk) = (1 + tZ ° + ^ + -) - + + + -) 

OO OO 

= £( E 4°'"^) f E E fr ( M »)'"' ( 54 ) 

n=0 k +... + k k =n n=0 

where the linear map from k + 1 dimensional vector space V to V is given by 
matrix MeGc GL(V) whose eigenvalues are zo, Zk- Using this observation 
several conclusions can be drawn. First, it should be clear that 

E *o° ■■■***= E x m = ir(M„). (55) 

k a + ...+k k =n m£riAnZ fc + 1 

Second, following Stanley [3, 27] we would like to consider the algebra of invari- 
ants of G. To this purpose we introduce a basis x ={ xq, Xk} of V and the 
polynomial ring R = C[xo, Xk] so that if / £ R , then M/(x) = /(Mx). The 
algebra of invariant polynomials R G can be defined now as 

R G = {f£R: M/(x) = /(Mx) - /(x) VM G G}. 

These invariant polynomials can be explicitly constructed as averages over the 
group G according to prescription: 

l G l MGG 

with \G\ being the cardinality of G. Suppose now that / € R a , then, evidently, 
/ e -R G = Avcf so that Av^f = Avcf = f ■ Hence, the operator Avg is 
indepotcnt. Because of this, its eigenvalues can be only 1 and 0. From here it 
follows that 

dim/f = ^ HM n ). (57) 



21 



Thus far our analysis was completely general. To obtain Eq.(9) we have to put 
Zo = ... = Zk = 1 in Eq.(8). This time, however we can use the obtained results 
in order to write the following expansion for the Poincare' polynomial [3, 27] 
which for the appropriately chosen G is equivalent to Eq.(lO): 

P(R G ,t) = E|^ E ^M n )r = £dim/„ G r. (58) 

n=0 ' ' MeG n=0 

Evidently the Ehrhart polynomial ^(k, n) = dim f£- To figure out the group G 
in the present case is easy since, actually, the group is trivial: G = l.This is so 
because the eigenvalues z , Zk of the matrix M all are equal to 1. It should 
be clear, however, that for some appropriately chosen group G expansion (19) 
is also the Poincare polynomial (for the Cohen -Macaulay polynomial algebra 
[3,27]). This fact provides independent (of Refs. [10,11]) evidence that both the 
Veneziano and Veneziano-like amplitudes are of topological origin. 



5.2 Connections with intersection theory 

We would like to strengthen this observation now. To this purpose, in view of 
Eq.(35), and taking into account that for the symplectic 2-form il = X^i=i dxi A 
dyi the n-th power is given byf2" = f2Af2A---AO =dx\ A dy\ A • • -dx n A dy n , 
it is convenient to introduce the differential form 

expfi = 1 + Q + ^yQ A Q + ifi A Q A ft H . (59) 

By design, the expansion in the r.h.s. will have only k terms. The form fi 
is closed, i.e. d£l = (the Liouvolle theorem), but not exact. In view of the 
expansion, Eq.(59), the D-H integral, Eq.(51) can be rewritten as 

I(k,f) = J exp(O), (60) 

kA 

where, following Atiyah and Bott [28], we have introduced the form Cl = O— f • x. 
Doing so requires us to replace the exterior derivative d acting on Q by d = 
d + (where the operator i(x) reduces the degree of the form by one) with 
respect to which the form is equivariantly closed, i.e.dft = 0. More explicitly 
we have ctt7 = dO, + i(x)Q — f • dx =0. Since dft — 0, we obtain the equation 
for the moment map : «(x)0 — f • dx =0 [51,58]. If use of the operator d on 
differential forms leads to the notion of cohomology, use of the operator d leads 
to the notion of equivariant cohomology. Although details can be found in 
the paper by Atiyah and Bott [28], more relaxed pedagogical exposition can 
be found in the monograph by Guillemin and Sternberg [60]. To make further 
progress, we would like to rewrite the two- form Q in complex notations [51]. To 
this purpose, we introduce Zj = pj + iqj and its complex conjugate. In terms of 
these variables ft acquires the following form : O = | dzihdzi. Next, recall 

[61] that for any Kahlcr manifold the fundamental 2-form ft can be written as 
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n = j Yljj hij(z)dzi/\dz,j provided that hij(z) = Sij +0(|z| 2 ). This means that 
in fact all Kahlcr manifolds are symplectic [26,54]. On such Kahler manifolds 
one can introduce the Chcrn cutrvature 2-form which (up to a constant ) should 
look like £1. It should belong to the first Chern class [19]. This means that, at 
least formally, consistency reqiures us to identify x^s entering the product f • x 
in the form f2 with the first Chern classes Cj, i.e. f • x = J2i=i fi°i- This fact was 
proven rigorously in the above mentioned paper by Atiyah and Bott [28] . Since 
in the Introduction we already mentioned that the Veneziano amplitudes can 
be formally associated with the period integrals for the Fermat (hyper)surfaces 
T and since such integrals can be interpreted as intersection numbers between 
the cycles on T [13,28,61] (see also Ref.[58], p. 72) one can formally rewrite the 
precursor to the Veneziano amplitude [10] as 

I = (t=|) ° • • • d J exp(^) | /i=0 vi= / dK(co) ro ■ ■ ■ (c d ) r « (61) 

'a a 

provided that ro + • • • + r^ = n in in view of Eq.(13). Analytical continuation 
of such an integral (as in the case of usual beta function) then will produce 
the Veneziano amplitudes. In such a language, calculation of the Veneziano 
amplitudes using generating function, Eq.(60), mathematically becomes almost 
equivalent to calculations of averages in the Witten-Kontsevich model [3 1-33] 9 . 
In addition, as was also noticed by Atiyah and Bott [28] , the replacement of the 
exterior derivative d by d = d + i(x) was inspired by earlier work by Witten on 
supersymmetric formulation of quantum mechanics and Morse theory [29] . Such 
an observation along with results of Ref . [60] allows us to develop calculations of 
the Veneziano amplitudes using supersymmetric formalism. 

5.3 Supersymmetry and the Lefshetz isomomorphism 

We begin with the following observations. Let X be the complex Hermitian 
manifold and let £ p+q (X) denote the complex -valued differential forms (sec- 
tions) of type (p,q),p + q = r, living on X. The Hodge decomposition insures 
that £ r (X)=J2p+ q =r £ p+q {X). The Dolbeault operators d and 8 act on S p+q {X) 
according to the rule d : £P+i{X) £P +1 ^(X) and 5 : £p+*(X) £P'i +1 (X) 
, so that the exterior derivative operator is defined as d = d + d. Let now 
<p p ,il> p £ P. By analogy with traditional quantum mechanics we define (using 
Dirac's notations) the inner product 

< <p p | V P >= J V P A *^ p , (62) 

M 

where the bar means the complex conjugation and the star * means the usual 
Hodge conjugation. Use of such a product is motivated by the fact that the 
period integrals, e.g. those for the Veneziano-like amplitudes, and, hence, those 

9 This fact is explained in more details in Ref. [34]. 
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given by Eq.(61), are expressible through such inner products [61]. Fortunately, 
such a product possesses properties typical for the finite dimensional quantum 
mechanical Hilbert spaces. In particular, 

< </> P I V> 9 >= CS P ,q and < *P P I <P P » 0, (63) 

where C is some positive constant. With respect to such defined scalar product it 
is possible to define all conjugate operators, e.g. d*, etc. and, most importantly, 
the Laplacians 

A = dd* + d*d, 

D = dd* + 8*8, (64) 
□ = 33* + 5*8. 

All this was known to mathematicians before Witten's work, Ref.[29]. The 
unexpected twist occurred when Witten suggested to extend the notion of the 
exterior derivative d. Within the de Rham picture (valid for both real and 
complex manifolds) let M be a compact Riemannian manifold and K be the 
Killing vector field which is just one of the generators of isometry of M, then 
Witten suggested to replace the exterior derivative operator d by the extended 
operator 

d s = d + si{K) (65) 

briefly discussed earlier in the context of the equivariant cohomology. Here s 
is real nonzero parameter conveniently chosen. Witten argues that one can 
construct the Laplacian (the Hamiltonian in his formulation) A by replacing A 
by A s = d s d* s + d* s d s . This is possible if and only if d 2 s = d* s 2 = or, since 
d 2 = sC(K) , where C(K) is the Lie derivative along the field K, if the Lie 
derivative acting on the corresponding differential form vanishes. The details 
are beautifully explained in the much earlier paper by Frankel, Ref.[62]. Atiyah 
and Bott observed that the auxiliary parameter s can be identified with earlier 
introduced f. This observation provides the link between the D-H formalism 
discussed earlier and Witten's supersymmetric quantum mechanics. 

Looking at Eq.s (64) and following Ref.s[14,51,58] we consider the (Dirac) 
operator 8 = 3+3* and its adjoint with respect to scalar product, Eq.(62). Then 
use of above references suggests that the dimension Q of the quatum Hilbert 
space associated with the reduced phase space of the D-H integral considered 
earlier is given by 

Q = kcr 8 — co kcr 8* . (66) 

Such a definition was also used by Vergne[14]. In view of the results of the 
previous section, and, in accord with Ref. [14], we make an identification: Q = 

We would like to arrive at this result using different set of arguments. To this 
purpose we notice first that according to Theorem 4.7. by Wells, Ref. [61], we 
have A = 2D = 20 with respect to the Kahler metric on X. Next, according to 
the Corollary 4.11. of the same reference A commutes with d, d*,8,8*,3 and 3*. 
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From these facts it follows immediately that if we, in accord with Witten, choose 
A as our Hamiltonian, then the supercharges can be selected as Q + = d+d* and 
Q~ = i (d — d*) . Evidently, this is not the only choice as Witten also indicates. 
If the Hamiltonian H is acting in finite dimensional Hilbert space one may 
require axiomatically that : a) there is a vacuum state (or states) | a > such 
that H| a >= (i.e. this state is the harmonic differential form) and Q + \ a >= 
Q- \a>=0 . This implies, of course, that [H,Q+] = [H,Q~] = 0. Finally, once 
again, following Witten, we may require that (Q + ) 2 = (Q~) 2 =H. Then, the 
equivariant extension, Eq.(65), leads to (Qf ) 2 = H+2is£(K). Fortunately, the 
above supersymmetry algebra can be extended. As it is mentioned in Ref. [61], 
there are operators acting on differential forms living on Kahlcr (or Hodge) 
manifolds whose commutators are isomorphic to sl 2 (C) Lie algebra. It is known 
[63] that all scmisimple Lie algebras are made of copies of sl 2 (C). Now we can 
exploit these observations using the Lefschetz isomorphism theorem whose exact 
formulation is given as Theorem 3.12 in the book by Wells, Ref.[61]. We are 
only using some parts of this theorem in our work. 

In particular, using notations of this reference we introduce the operator L 
commuting with A and its adjoint L* = A . It can be shown, Ref. [61], p. 159, 
that L* = w * L* where, as before, * denotes the Hodge star operator and the 
operator w can be formally defined through the relation ** = w, Ref.[61] p. 156. 
From these definitions it should be clear that L* also commutes with A on the 
space of harmonic differential forms (in accord with p. 195 of [61]). As part of 
the preparation for proving of the Lefschetz isomorphism theorem, it can be 
shown [61], that 

[A, L] = B and [B, A] = 2A, [B, L] = -2L. (67) 

At the same time, the Jacobson-Morozov theorem, Ref. [36], and results of 
Ref. [63], p. 37, essentially guarantee that any sl 2 (C) Lie algebra can be brought 
into form 

[h a ,e a ] = 2e a , [h a , f a ] = -2f a , [e a ,f a ] = h a (68) 

upon appropriate rescaling. The index a counts thenumber of sh^C) algebras 
in a scmisimple Lie algebra. Comparison between the above two expressions 
leads to the Lie algebra endomorphism, i.e. the operators h a , f a and e a act 
on the vector space {v} to be described below while the operators A, L and B 
obeying the same commutation relations act on the space of differential forms. 
It is possible to bring Eq.s(67) and (68) to even closer correspondence. To this 
purpose, following Dixmier [64], Ch-r 8, we introduce operators h — ^ Q a a h a , 
e = J2 a ^a e ai f = J2 a c afa- Then, provided that the constants are subject to 
constraint: b a c a = a a , the commutation relations between the operators h, e 
and / are exactly the same as for B, A and L respectively. To avoid unnecessary 
complications, we choose a a = b a = c a = 1. 

Next, following Serre, Ref. [35], Ch-r 4, we need to introduce the notion of the 
primitive vector (or element). This is the vector v such that hv=Xv but ev = 0. 
The number A is the weight of the module V x = {v 6 V hv=Xv}. If the 
vector space is finite dimensional, then V — J2\V X - Moreover, only if V x is 
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finite dimensional it is straightforward to prove that the primitive element does 
exist. The proof is based on the observation that if x is the eigenvector of h 
with weight A, then ex is also the eigenvector of h with eigenvalue A — 2, etc. 
Moreover, from the book by Kac [65], Chr.3, it follows that if A is the weight of 
V, then A— < A, > on is also the weight with the same multiplicity, provided 
that < A,o^ >G Z. Kac therefore introduces another module: U = J2kez 
V x+kai . Such a module is finite for finite Weyl-Coxeter reflection groups and 
is infinite for the affine reflection groups associated with the affine Kac-Moody 
Lie algebras. 

We would like to argue that for our purposes it is sufficient to use only 
finite reflection (or pseudo-reflection) groups. It should be clear, however, from 
reading the book by Kac that the infinite dimensional version of the module U 
leads straightforwardly to all known string-theoretic results. In the case of CFT 
this is essential, but for calculation of the Veneziano-like amplitudes this is not 
essential as we are about to demonstrate. Indeed, by accepting the traditional 
option we loose at once our connections with the Lefschetz isomorphism theorem 
( relying heavily on the existence of primitive elements) and with the Hodge 
theory in its standard form on which our arguments are based. The infinite 
dimensional extensions of the Hodgc-dc Rham theory involving loop groups, etc. 
relevant for CFT can be found in Ref. [66] . Fortunately, they are not needed for 
our calculations. Hence, below we work only with the finite dimensional spaces. 

In particular, let v be a primitive element of weight A then, following Serre, 
we let v n = ^e n v for n > and V-\ — 0, so that 

hv n = (A - 2n)v n 
ev n = (n + l)v n+ i 
fv n = (A - n + l)u n _i 

Clearly, the operators e and / are the creation and the annihilation operators 
according to the existing in physics terminology while the vector v can be in- 
terpreted as the vacuum state vector. The question arises: how this vector is 
related to the earlier introduced vector | a >? Before providing an answer to 
this question we need, following Serre, to settle the related issue. In particular, 
we can either: a) assume that for all n > the first of Eq.s(69) has solutions 
and all vectors v,v\ 1 V2 , are linearly independent or b) beginning from 
some m + 1 > 0, all vectors v n are zero, i.e. v m ^ but v m+ \ = 0. The first 
option leads to the infinite dimensional representations associated with Kac- 
Moody affine algebras just mentioned. The second option leads to the finite 
dimensional representations and to the requirement A = m with m being an 
integer. Following Serre, this observation can be exploited further thus leading 
us to crucial physical identifications. Serre observes that with respect to n = 
Eq.s(69) possess a ("super") symmetry. That is the linear mappings 

e m :V m ^ V- m and f rn : V~ m -» V m (70) 

are isomorphisms and the dimensionality of V m and V~ m are the same. Serre 
provides an operator (the analog of Witten's F operator) 9 = exp(/) exp(e) exp(— /) 



(69) 



26 



such that 6 ■ f = —e ■ 6, 6 ■ e = —9 ■ f and 9 ■ h = —h ■ 9. In view of such an 
operator, it is convenient to redefine h operator : h — > h = h — A. Then, for 
such redefined operator the vacuum state is just v. Since both L and L* = A 
commute with the supersymmetric Hamiltonian H and, because of the group 
endomorphism, we conclude that the vacuum state | a > for H corresponds to 
the primitive state vector v. 

Now we are ready to apply yet another isomorphism following Ginzburg [36] , 
Ch-r. 4, pp 205-206 10 . To this purpose we make the following identification 

d d ( d d \ 

e * u+1 du > fi ^ u du^> hi ^ 2 { u+1 duTi ~ u du)> (71) 

i = 0, m. Such operators are acting on the vector space made of monomials 
of the type 

Vn - ?n = f C ■ • • t n k " , (72) 

n !ni! • • • n k \ 

where no + ... + = n . This result is useful to compare with Eq.(61). 
Eq.s (69) have now their analogs 

hi * T n {i) = 2(n i+ i - rii)T n {i), 

e l *T n {i)^2n l T n {i + l), (73) 
fi * F n {i) = 2n i+ iF n (i - 1), 

where, clearly, one should make the following consistent identifications: m(i) — 
2n(i) — 2 (rii+i —rii) , 2rii — n(i) + 1 and m(i) — n(i) + 1 = 2rii + i. Next, 

we define the total Hamiltonian: h = J2i=o^i Uso tnat Si=o m (0 = n ' an( i 
then consider its action on one of the wave functions of the type given by 
Eq.(72). Since the operators defined by Eq.s(71) by design preserve the total 
degree of monomials of the type given by Eq.(72) (that is they preserve the 
Veneziano energy-momentum codition) , we obtain the ground state degeneracy 
equal to *P(fc, n) in agreement with Vergne, Ref.[14], where it was obtained using 
different methods. Clearly, the factor *P(fc,n) is just the number of solutions in 
nonnegative integers to n + ... + n k = n, Rcf.[53], p. 252. 



10 Unfortunately, the original sourse contains absolutely minor mistakes. These are easily 
correctable. The corrected results are given in the text. 

n The physical meaning of h is discussed in some detail in our earlier work, Ref.[ll]. 
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